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I. INTRODUCTION 



A classical gas of charged point particles is non- magnetic, by Van Leeuwen's Theorem [|l|. 
But quantum mechanically, free spinless Bose particles in a uniform magnetic field show 
diamagnetism ||^. One finds that the energy of the system in a magnetic field is higher 
than in the absence of the field. Simon proved quite generally that the ground state of 
non-relativistic spinless Bosons interacting through an arbitrary potential always increases 
in a magnetic field. He went on to extend this result by showing that the free energy in 
the presence of a magnetic field is always greater than the free energy in the absence of a 
magnetic field at all temperatures An alternative proof of this result is given in [^]. All 
this work which deals with systems with a finite number of degrees of freedom suggests that 
diamagnetism is a universal property of Spinless Bosons. In field theory (which describes 
systems with an infinite number of degrees of freedom) charged spinless Bosons are described 
by complex scalar fields. One might therefore expect that charged scalar fields would also 
show diamagnetic behaviour. With this motivation we study the magnetic behavior of scalar 
field theories. 

The paper is organised as follows. The first part deals with finite temperature free scalar 
field theory in the presence of an external homogeneous magnetic field. Here we explicitly 
calculate the partition function and the free energy as a function of the applied magnetic 
field. This expression is formally divergent. Using a suitable regularization scheme we 
compute the difference in the free energy (with and without the magnetic field) and obtain 
a finite answer. This difference is also shown to be positive, thus establishing the diamagnetic 
behaviour of free charged scalar fields. 

We then move on to interacting scalar field theory in the second part. Here, we cannot 
evaluate the partition function explicitly. Nevertheless we prove the universal diamagnetism 
of scalar fields by assuming a finite momentum cutoff in the theory. If the theory is renor- 
malizable, then one can take this cutoff to infinity while maintaining finiteness of all physical 
quantities. In both cases the results obtained are exact. 
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II. FREE CASE 



In this section we calculate the free energy of free scalar fields in the presence of an 
external uniform magnetic field. For ease of presentation we work in two spatial dimensions. 
The interesting physics takes place in the plane normal to the applied field. Generalization 
to higher dimensions is straightforward. 

Let $ be a complex scalar field which describes charged spinless Bosons. The Lagrangian 
density of a free charged scalar field in the presence of a constant homogeneous external 
magnetic field is given by 



L = {p^^)*{p^^) - m'($*$) (1) 



where = 0, 1, 2, 



D^^d^- leA^ (2) 



and m and e are the mass and charge respectively. (We set ^ = 1 and c — 1). Now, we 
write the complex field in term of two real fields $i and $2- 

$ = 1= — = 1= — (3) 

V2 \/2 ^ ' 

This theory has a global U(l) symmetry and therefore a conserved Noether charge Q, given 

by 

g = y d2x(ni$2-$in2) (4) 

where 

= do^i (5) 

The Hamiltonian density of the system is given by 

^ = ^(n? + n^) + ^(v$if + ^(v$2f + 

\{m' + e'A'){^l + ^l)-j.A, (6) 



where the current density j is given by 

j = _e($iV$2 - $2V$i). (7) 

We now suppose that the external magnetic field is uniform in the x — y plane. 
We choose the temporal gauge { Aq = 0). The constant magnetic field B is 

B=da:Ay-dy (8) 

where A(x) is independent of t. 
The action of this theory is 

S = j d^xdTC{^,^\A), (9) 

where r is the imaginary time variable which runs from to /3 (=1/(A;bT)), the inverse 
temperature. The action defined above is quadratic and so the partition function can be 
evaluated exactly. As is usual in finite temperature field theory 0, we impose periodic 
boundary conditions for Bosonic fields 

$(x,0) = <l>(x,/?). (10) 

Now, the partition function of this theory can be written as 



z{B) = j i?[ni]i?[n2] j P[$i]i)[<i>2]exp 



drd X 



iHi^ + in2^ - <^'2, Hi, Hs) + /i($2ni - $in2) 



(11) 



Here [i is the chemical potential associated with the conserved charge. We pick the gauge 
in which the vector potential A is (-By, 0), and expand the complex scalar field in terms 
of modes adapted to the present situation. These modes solve the Klein-Gordon equation 
in an external magnetic field. The eigenfunctions are labelled by one discrete (/) and one 
continuous 'Px quantum number and the spectrum depends on I only. In the gauge we 
choose, the modes are plane waves in the x direction and harmonic oscillator (i.e. gaussian) 
wavefunctions in the y direction. 



The spectrum is given by 

u;2 = + (2/ + l)e5, 1 = 0,1,2 oo (12) 

The degeneracy of these states for fixed / is eAB /27r , where A is the area of the system. 
So, these modes can be thought of as quantized harmonic oscillators. Expanding the fields 
$1 and $2 in these modes the system reduces to a collection of harmonic oscillators with 
frequency ui. 

By standard manipulations H we get the free energy as 



F{B) = - In Z(B)//3 = 27reA5 ^ 

1=0 L 



u;/ + ln(l - exp{-P{uJi - fi)) 



(13) 



The first term in the square brackets corresponds to the zero point fluctuation of the 
vacuum and the other two terms are finite temperature contributions. We will first compare 
the free energy of the system with and without the magnetic field at zero temperature. 

The free energy of the system in presence of the magnetic field at zero temperature is 
given by 

oo 

Fo{B) = 27r A eBY. coi, (14) 

1=0 

where uf = rn? + {21 + l)eB. Obviously, this sum diverges. In order to obtain a finite answer 
we need to impose a cutoff L in the sum (|14|). Then the free energy becomes 

L 

Fo{B, L) = 2nA eB ^ ui. (15) 

1=0 

The free energy in the absence of the magnetic field is also given by a divergent expression 

Fo(0) = '271 A pdp ^{pI+pD + tv? (16) 
We regularize this expression by imposing a cutoff A. Thus the free energy ([l6|) becomes 

Fo(0, A) = 271 A £ pdp ^{pI + pD + m2 (17) 



In order to compare the free energies in equations (|T5D and (|T^, we choose the cutoffs L and 
A in such a way that both systems have the same number of modes. Such a procedure can 
be justified on physical grounds if one imagines turning on the magnetic field adiabatically. 
Counting the modes upto the L-th Landau level we find 

L 

2^ A eBY, 1 = 2vrA eB{L + 1) (18) 

1=0 

Similarly, for the momentum cutoff upto A we get the modes without the magnetic field as 

/•A 

27rA / pdp = 7r AA'^ (19) 
Jo 

Equating these gives us 

A^ = 2 eS (L + 1) (20) 

Now, the free energy in absence of the magnetic field depends on magnetic field through the 
momentum cutoff and is given by 

„A(ij) 

Fo(0, L) = 27tA / pdp ^p^ + m? (21) 

•J 

The difference between the two free energies is given by 

AF(5, L) = FoiB, L) - Fo(0, L) (22) 

We define G{B) = Fo{B,L)/27rA, H{B) = Fo(0,L)/27rA and K{B) = G{B) - H{B). 
Numerically evaluating these sums for large but finite L and plotting them in figure 1 and 
figure 2 shows that K{B) is the difference between two large quantities. As the cutoff L 
goes to infinity then K{B) becomes the difference between two infinities. In this limit we 
find that K{B) tends to a finite value. Thus, the susceptibility at zero temperature in the 
relativistic case is non-zero. This vacuum susceptibility can be interpreted as due to virtual 
currents. 

We now show analytically that AF{B) is positive i.e. the vacuum is diamagnetic. Note 
that 
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AF{B) = F{B)-F{0) = J2 ai{B,m) 

1=0 



(23) 



where ai{B,m) is given by 



ai{B, m) = eB 



m2 + {21 + l)eB - da ^m? + 2(/ + a)eB 



(24) 



Introducing a dimensionless quantity p = ^ the above equation becomes 



+ {21 + l)p - rfa \/l + 2{l + a)p 



(25) 



The positivity of ai{p) for each / can be proved geometrically. Defining zi = {1 + 2lp) /2p 
and f{a) = y/zi + a, the coefficient ai{p) can be rewritten in terms of q(p) as 



(26) 



Since, the function f{a) is convex, the area under the tangent drawn at a = 1/2 is greater 
than the area under the curve (see figure 3). This shows that ci{p) is positive. To show the 
convergence of the sum (^) we note that 



/(1/2) - 



/(0) + /(l) 



/ da f{a) < 
Jo 

Now, applying mean value theorem twice one can easily show that 

/' da f{a) < - ^ 3/2 
JO 10(2;^ + a) 



(27) 



(28) 



Thus the coefficient q(p) is positive for each / and the sum converges, hence the diamagnetic 
inequahty is established. 

Massless Limit: The magnetisation at zero temperature in the zero mass limit is given 



by 



M(B) 



(29) 



So, the susceptibility in this zero mass limit is given by 

X{B) ^ 



(30) 



which diverges [^Q as B goes to zero. This feature of the susceptibihty has aheady been 
noticed in the magnetised pair Bose gas p]. So, in this case the external field will be totally 
expelled. This happens because of large number of virtual particle and antiparticle produced 
in the ground state so that the overall diamagnetism of the system is high enough to totally 
expel the external field. 

Another interesting point is that there exists a critical magnetic field below which mag- 
netic field will be totally expelled. This critical field can be estimated as follows. The 
effective magnetic field can be defined as -Be// = B + 271 M Here, the magnetisation M 
varies as — const\fB. Therefore, there exists a critical magnetic field Be where the effective 
field Beff vanishes. 

Now, for the finite temperature case one can regulate the free energy through the same 
mode matching regularisation method. Finally, one can write down the free energy difference 
in dimensionless form as before 

oo 

AF(5) = F{B) - F(0) = blip, C) (31) 



where, 



1 



g{p,l, 1/2) - / da g{p,l,a) 







(32) 



The dimensionless variables are defined as 5 = jSm and ( = jSp. 
The coefficient g{p,l, a) is given by 



g{p, I, a) = log[l- exp{-5{Jl + 2{l + a)p - 0) + 



log [I - exp{-6{^l + 2{l + a)p + C)) j • (33) 

Now, defining zi = ^ ^^2^^'^^ we can write the equation (^) 

g{p, /, a) = log (l - exp{-{^/zi + a - C))) + 

log (1 - exp{-{^zi + a + Q)) ■ (34) 

The function g{p, /, a) is convex and so the zero temperature argument applies unchanged. 
It follows that the free energy satisfies the following inequality 
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F(B) > F(0) 

Thus the response of the system to the magnetic field will be diamagnetic. 



(35) 



III. INTERACTING CASE 

In this section we want to extend the diamagnetic inequality to the self-interacting field 
theory case including the dynamical interaction between scalar fields. The partition function 
of this charged self-interacting field theory in the presence of the magnetic field can be written 

as 

Z(B) = J y'p[$]P[$*]exp(-5($,<l>*,A)), (36) 

where the action S is defined as 

d^x dr [{D^^){D^'^)* + m2($*$) + $*)] . (37) 

The action is not quadratic and Z{B) cannot be evaluated in closed form. Nevertheless, 
we show that the response of the system to an external magnetic field is diamagnetic. Since 
the formal expression for the partition function may not exist (the integrals may not exist) 



we impose a cut off in momentum space. The functional integral in (|36|) signifies that one 
only integrates over those field configurations whose Fourier transforms have support within 
a sphere of radius A in momentum space. The partition function then explicitly depends on 
A. We do not explicitly indicate the A and /i dependence of Z{B, A,/i) below. 

We divide the action into two parts Sq and Sint, where Sq is the action in the absence of 
the external field. 

S = So + Sint, (38) 

where 

So= / f d^xdT\{df,<^){d^<l>y +m\^*^)+V{<^*,<^)], (39) 
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Jint 



d'^xdr 



-ie(5/,$)(A^<l>*) + ie(A^$)(a'^$*) + e2(A- A)(<l><l>*) . (40) 



Notice that exp{—So) is a positive measure on the space of field configurations. The ra- 
tio Z{B)/Z{0) can therefore be regarded as the expectation value of exp{—Sint)- Since 
exp{—Sint) is an oscillatory function whose modulus is less than or equal to 1, we conclude 
that 

^^^^ exp-5i„4 >< 1 (41) 



Z(0) 



This implies that 



F{B) > F(0) (42) 



In this derivation, we have not assumed any form for the vector potential. So, the 
result derived above is true for both homogeneous or inhomogeneous magnetic fields of any 
strength. Since (3 is arbitrary, the result holds at all temperatures. The argument presented 
here works for any arbitrary interaction V^( $*$).( Generally, it is assumed that V($*$) is a 
smooth function, for instance, a polynomial [Q). 

Upto now we have considered the cases of charged scalar fields interacting through a 
potential. It is also possible to consider interaction mediated by a djTiamical electromagnetic 
field A^. The fields in the system are now $ (charged scalar fields) and A^. If one applies 
an external magnetic field Ag^t then the full Lagrangian is given by 

C = -^F^ + {D^<^)*{D^'<^) - m^{<t>*<^) - V{<^*, $) (43) 

where D^ = d^- ieA^^* - ieA^ and F^,^ = d^A^, - d^A^,. 

The argument given above can be modified as follows. The definition of 5*0 changes 
slightly while Sint remains the same. 

1 



d xdr 



4 



F' + {d^ - ^eAf,^y{^^' + ieA^"^) + m\^*^) + 1/($*<l>) 



(44) 



and 
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Jint 



d'^xdr 



(45) 



Again one can repeat the same argument by noting that exp{—So) is a positive measure 
and the ratio Z{B)/Z{0) as an expectation value of exp{—Sint) to estabhsh the diamagnetic 
inequahty. 



IV. CONCLUSIONS AND PERSPECTIVES 

The response of a system to an electric field is completely different from its response to 
a magnetic field. The basic difference between the responses of a system on application of 
an electric field or a magnetic field lies in the Hamiltonian of the system. 

The Lagrangian of a system in the presence of an electric field can be written as 

£ = (Do$)*(£'o'^) - (V$)*(V$) - m2($*<l>) - \/(<l>*$) (46) 

where 

Do = do- ieAo (47) 

For statistical mechanics to make sense, the Hamiltonian H must be independent of time. 
Then it follows that 

n = (n*)(n) + (v<i>)*(v$) + m2($*$) + r($*<i>) - le [(n*)(Ao$) - (n)(Ao$*)] (48) 

The electric field appears in the Hamiltonian through the linear vector potential Aq term. 
Now, from finite temperature second order perturbation [|l3,|ll|] theory, one can show easily 



that the free energy of the system always decreases with the electric field. Hence, the 
dielectric susceptibility is always positive in thermal equilibrium. 

But in the case of a magnetic field the Hamiltonian contains both linear and quadratic 
terms in A. The net effect of an applied magnetic field is not a priori clear. However, as our 
analysis makes clear, for charged scalar field theories the net effect is always diamagnetic. 
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In case of Spinless Bosons, there is no Zeeman term coupled with a magnetic field 
and hence the system consisting of Spinless Bosons always has higher energy in a magnetic 
field than without the magnetic field. It has been already pointed [^,^|13|,|14| out that there 
is no corresponding theorem for Fermions. 

Let us consider some illustrative examples of Spinless Bose systems. One obvious example 
in the laboratory is the Cooper pair formation in superconductors which shows perfect 
diamagnetism (known as the Meissner |jl5| effect) below the critical temperature. Cooper 



pairs also exist in Neutron Stars ||T6[ where the magnetic field is very high compared to any 
laboratory field. Of course, the operators which create and destroy Cooper pairs are not 
strictly Bose operators, so this is only an analogy. It has been shown explicitly by J. Daicic 
et. al. that the magnetised pair Bose systems are relativistic superconductors. These 
systems are not covered by previous analyses which apply to non-relativistic quantum 

mechanical systems. Pions would be a suitable candidate for apphcation of this theory with 
7r+ and tt^ regarded as the particles and antiparticles. The choice is also motivated by the 
fact that these pions are massive (mc^ = 139.5673 Mev), obey Bose-Einstein statistics and 
that they possess no spin. It is also well known that hard core Bosons in any dimension 
on any lattice show a preference for zero fiux. 

Throughout this paper we have considered two spatial dimensions both for the free and 
interacting cases, but the formalism easily generalizes to higher dimensions. In summary, 
we have shown exactly, that charged scalar fields at all temperatures are diamagnetic. 
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FIGURES 

FIG. 1. The functions G{B) and H{B) have been plotted against B for L = 3 and m = .1. For 

simpUcity we have drawn the figure for B instead of eB. 

FIG. 2. The function K{B) = G{B) — H{B) has been plotted against B. The values of L and 
m are the same as figure 1. Note that the scale in this figure in y direction is much expanded 
compared to figure 1. 

FIG. 3. The full line is the curve f{a) = \/ .1 + a, the dashed line is the tangent to the above 
curve at a = .5. It is straightforward to see that the area under the curve is less than the area 
under the tangent. Also the area under the tangent is same as the area of the rectangle defined by 
the dotted line. The area of this rectangle is given by /(1/2) x 1 = /(1/2). Hence, the positivity 
of a(.l) is proved. This can be generalised to any positive value of zi. 
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